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ABSTRACT 

We present observations, analysis and results for the first-year operation of AMiBA, an interferomet- 
ric experiment designed to study cosmology via the measurement of Cosmic Microwave Background 
(CMB). AMiBA is the first CMB interferometer operating at 3 mm to have reported successful re- 
sults, currently with seven close-packed antennas of 60-cm diameter giving a synthesized resolution of 
around 6^ During 2007 AMiBA detected the Sunyaev-Zel'dovich effects (SZE) of six galaxy clusters 
at redshift 0.091 < z < 0.322. An observing strategy with on-off-source switching is used to minimize 
the effects from electronic offset and ground pickup. Planets were used to test the observational ca- 
pability of AMiBA and to calibrate the conversion from correlator time-lag data to visibilities. The 
detailed formalism for data analysis is given. We summarize our early tests including observations of 
planets and quasars, and present images, visibility profiles, the estimated central coordinates, sizes, 
and SZE amplitudes of the galaxy clusters. Scientific implications are summarized. We also discuss 
possible systematic effects in the results. 

Subject headings: cosmic microwave background — cosmology: observations — methods: data anal- 
ysis — galaxies: clusters 



L INTRODUCTION 

The Cosmic Microwave Background (CMB) has been 
used as a window through which to study not only 
the early Universe but also its evolutionary history. 
The challenge of such observations arises from the fact 
that the CMB carries information spanning a period 
of about 14 billion years, so that specially-designed in- 
struments and analysis methods are required to sepa- 
rate the effects originated from different physical pro- 
cesses. One recent focus is the measurement of the 
Sunyaev-Zel'dovich effect s (SZE) result ings; from the hot 
g:as in galaxy clusters dSunvaev and ZePdovich I 1 19721 : 
iBirkinshaw fc Lancaster! l2007l ). Several CMB exper- 
i ments are dedi cated t o this purpose, su ch as ACT 

j^osowskv'2003'), AMI (jKneissl et al\m^ , APEX-SZ 
ialverson et al 2m§t)^ OVRO /BIMA-SZE (J ^aRoQue 
306), SPT (Ruhl 2004). SZA (iMuchovei et al\\2mi^ . 
and OCRA-p (Lancaster et "ZI I2007D . Here we report 
the first results of SZE observations with the Y. T. Lee 
Array for Micro wave Background Anisotropy (AMiBA; 
iHo et a/.l[2008ah . 

AMiBA is devoted to CMB observation with par- 
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ticular emphasis on detecting SZE clu sters, the CMB 
anisot ropy, and poss i bly cosmic defec t s (iLin et a l! '2004'; 
Umetsu ei^ g/J l2004l : IWu et al\ l2004l : iHo et all 2008b; 
Wu et g/T 2008af )! The instrument is described in 
Ho et al\ (|2008af ). IChen et all (|2008[ ). and iKoch et all 
(|2008a[ ). In our initial configuration with seven close- 
packed antennas of 60-cm diameter, we concentrated on 
pointed observations of the SZE and the measurement of 
the CMB temperature power spectrum (|Ho et g/.ll2008bl : 
IWu et g/.l[2008al 

At the AMiBA operating frequency (86-102 GHz; - 3- 
mm wavelength) the SZE is expected to induce a decre- 
ment in CMB intensity as compared with the undistorted 
Planckian spectrum (see Sec. [2j). This provides a pow- 
erful tool for discriminating between galaxy clusters and 
other astronomical sources, because the latter normally 
emit photons characteristic of much harder spectra than 
the CMB, and therefore always induce an increment in 
intensity rather than a decrement. By measuring the 
SZE intensity deficit and its profile, we hope to probe 
not only the cluster physics but also the related cosmic 
origins. 

This paper is organized as follows. In Section 2 we de- 
scribe the AMiBA telescope and its relation to the SZE. 
In Section 3, we summarize our initial observations in 
2007 with emphasis on the observing strategy and target 
selection. In Section 4, we describe our analysis formal- 
ism and procedures for obtaining the calibrated visibil- 
ities from the raw time-lag data, and the SZE images 
and profiles from the calibrated visibilities. The main 
results are also presented here. In Section 5 we discuss 
possible systematic errors. In Section 6, we summarize 
the scientific implications of the se results, includ ing the 
estimation of Hubble constant (|Koch et g/.ll2008b) . the 
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scaling rela tionship between t he SZE and X-ary derived 
properties (|Huang et al\ I2QQ8I) , and the baryonic frac- 
tion when our SZE data are jointly ana lyzed with the 
Subaru lensing data (|Umetsu et a/.l [20081 ). Brief conclu- 
sions are given in Section 7. We have other c ompanion 
pape rs investigating the system perform ance (|Lin et al\ 
I2008D . the foreground (iLiu et al\ I2008D . and data in- 
tegrity (|Nishioka et a/.ll2008[ ). 

2. AMIBA TELESCOPE AND SZE 

AMiBA is an interferometric experiment initiated in 
Taiwan in 2000 and dedicated on Mauna-Loa (3400 m 
in elevation), Big Island, Hawaii on October 3, 2006. It 
has dual-channel receivers operating at 86-102 GHz, de- 
signed to have full polarization capabilities (|Chen et al\ 
[2OO8). Currently it has sev en close-pac ked Cassegrain 
antennas of 60-cm diameter (|Koch et a/.l 12006) giving a 
synthesized resolution of 6 arcminutes (see Sec. 13. 2p . ex- 
pandable to a total of 19 elements with a synthesized 
resolution of about 2 arcminutes. The project has been 
funded for an expansion to 13 elements with dishes of 1.2- 
m diameter. The 13-element system is expected to start 
operating in the early 2009. Its cap ability in studying 
the SZE science is investigated by Mol nar et a l. (2008). 

The receiver-antenna elements are reconfigurable and 
co-mounted on a s ix- meter platform, which is driven by 
a hexapod mount (|Koch et al\\2008sl\ ). Each element has 
a cooled heterodyne receiver, consisting of HEMT am- 
plifiers of 46 dB in amplification, subharmonic mixers, 
and 2-18 GHz IF amplifiers. For each baseline, the sig- 
nals from two dual-channel receivers are cross-correlated 
in an analogue 4-lag correlator, whose time-lag outputs 
are convertible to two complex visibilities correspond- 
ing to the upper and lower frequency bands (see Sec.[4|). 
The cross correlation between the L and R polarization 
modes (the dual channels) of a pair of receivers enables 
the measurement of the four Stokes parameters, T, Q, 
t/, and V. Currently AMiBA operates with only two 
cross-polarization modes of LL and RR, focusing on the 
measurement of T. The typical recei ver noise te mpera- 
tures are between 80 and 110 K (L in et a/.ll2008f ). 

In 2007, the seven-element AMiBA focuses on targeted 
SZE observations. The theoretically expected SZE in 
fiux density is 

A/sze(/) = A/tszE(/,^) + A/kszE(/,r,^p), (1) 

where A/tszE and A/kSZE are the thermal and kinematic 
SZE respectively, defined as 

A/tszE(/, y) = loy [g{x) + 5,M. ^e)] , (2) 

A/kSZE(/,r,^p) = -/o/3rMx). (3) 

Here / is the observing frequency, Iq = 
2(kTcM^f / {hcf ^ 2.7 X 10^ Jy Sr-^ for a CMB tem- 
perature of TcMB = 2.725K, y = kaT / {rrieC^) J TqUqcII is 
the Compton parameter, g{x) = h{x)[x/ tanh(x/2) — 4] 
with X = hf /{kTcMB)^ ^rei is a relativistic correction, 
/3 = Vp/c is the peculiar radial velocity in units of 
speed of light, r = ax / riedl is the optical depth, 
h{x) = x^e^/{e^ — 1)^, c is the speed of light, ttIq is 
the electron mass, Uq is the electron number density 
in the cluster, Tq is the electron temperature, k is the 
Boltzmann constant, h is the Planck constant, and ax 
is the Thomson cross section. 



For the thermal effects, the A/tszE has a maximum 
decrement with respect to the CMB Planckian spectrum 
/cmb(/) = hx^/{e'' - I) 3i f ^ 100 GHz. This opti- 
mal frequency for SZE observations is well covered by 
the AMiBA operating frequency of / = [86, 102] GHz. 
For AMiBA g{f) ^ —3.4 at the center frequency / = 
94 GHz, with a variation of less than ±10% for the range 
of [86, 102] GHz. Therefore for a typical massive cluster 
with y - 10-"^, Te - 10 keV, r - ym^c^/{kT^) ^ 5 x 
10~^, Vp = 1000 km/s, and an angular size of = 6^ we 
expect the SZE signals to be A/tszE7r(6>/2)2 - -200 mJy 
and |A/kszE|7r(i9/2)2 - 20 mJy for AMiBA, indicating 
|A/kszE/A/tszE| ^ 10 %. The fractional contribution 
from the relativistic correction ^rei in such case is about 
7 %. Thus we expect the AMiBA SZE signals to be 
dominated by the thermal effects, with a feature of decre- 
ment in fiux density. We also note that |/cmb(94GHz)| ^ 
2.9 X 10^ Jy St-\ so that | A/tszE//cMB | ^ 3 x lO""^. 

The operating frequency of [86, 102] GHz also has the 
advantage that the SZE signals are less affected by the 
point sources when compared with the observations at a 
lower frequency, which is more typical for interferomet- 
ric ex periments such as A MI (15 GHz; iKneissl et al\ 
[200lh . SZA (27-35 GHz : iMuchovei et al\ [2007h. and 
OVRO/BIMA (30 GHz; iReese et a/.l 120021. This is be- 
cause typically the point sources have power-law spectra 
of flux density with negative spectral indices, so that 
they contribute less towards higher frequencies. We note 
that although SZA also has instruments at 90-98 GHz, 
their baselines are about ten times longer than ours pro- 
viding information at much smaller angular scales that 
is complementary to what AMiBA can supply. Another 
advantage of the AMiBA frequency range is that it also 
suppresses the Galactic synchrotron radiation and dust 
emission. Detailed investigation about the AMiBA fore- 
ground is presented in Liu et aQ (|2008l ). These features 
in frequency and resolution make AMiBA a unique CMB 
telescope, complementing the capability of other CMB 
projects. 

3. OBSERVATIONS 

3.1. Primary Beams 

Before the 60-cm antennas were mounted onto the 
platform, their beam patterns were individually mea- 
sured in the far field (~ 100 m) by scanning a wide- 
band white source with a system retrofitted from a 
commercial equatorial mount (Koch et al 2006). Fig- 
ure [TJa) shows an example. An azimuthal average of 
the beam pattern (Figure [Tfb)) reveals a Full-Width-at- 
Half-Maximum (FWHM) of about 23 arcminutes. The 
first side lobe is about 20 db below the primary peak 
and the beam profile i s as per designed. The index of 
asymmetry (jWu et a/.ll200ll ) is calculated showing that 
the discrepancy of the beam pattern from its azimuthal 
average is below 10% within the FWHM. All beam pat- 
terns of the antennas are well fitted by a Gaussian beam 
of 23 arcminutes within the 2-a region for < 10% error. 

3.2. u-v Coverage and Resolution 

In 2007, AMiBA was used with the seven 57.6-cm 
diameter antennas close-packed on the platform, pro- 
viding 21 baselines of three different lengths = 
60.6 cm, 105.0 cm, and 121.2 cm. Each base- 
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Fig. 1. — (a) Left: the measured primary beam of one of 
the 60-cm antennas, (b) Right: azimuthally averaged beam 
shape (sohd) and a Gaussian fit (dashed); the dotted hues 
indicate the FWHM. In both plots the intensity is in units of 
db normahzed to the peak value. 
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Fig. 2. — The u-v coverage of a typical observation by seven- 
element AMiBA (left) and its corresponding noise-weighted 
synthesized beam (the dirty beam; right). In the left plot, 
the color scale indicates the relative sensitivity in u-v space, 
while the dots mark the u-v modes sampled by the typical 
AMiBA observation. 

line operates at two frequency channels centered at 
/ = 90 and 98 GHz, each with a band width 
of 8 GHz. Thus the three baseline lengths corre- 
spond to six multipole bands £ = 2'k\Ju^ -|- = 
27rL^//c = [1092,1193], [1193,1295], [1891,2067], 
[2067,2243], [2184,2387], [^2387,2590]. Th ese baselines 
have a six-fold symmetry (|Ho et a/]l2008af ). To achieve 
better u-v coverage and thus a better imaging capabil- 
ity, we observe each target at eight different polariza- 
tion angles of the platform with an angular step of 7.5 
degrees, uniformly covering azimuthal angles in the u-v 
plane with 48 discrete samples for each of the three differ- 
ent baseline lengths at each of the two frequency bands. 
This makes a total of 288 discrete u-v samples. Figure [2] 
shows such typical u-v coverage, with the correspond- 
ing noise-weighted point spread function from a typical 
observation, the so-called 'dirty beam'. The FWHM of 
the dirty beam is about 6 arcminutes, defining the syn- 
thesized resolution of the seven-element AMiBA. When 
using an equal-weighted (rather than the current noise- 
weighted) scheme to construct the dirty beam, we still 
obtain a similar resolution of about 6 arcminutes because 
the noise level for each baseline is about the same. 

3.3. Observing Strategy 

To minimize, if not completely remove, the effects from 
ground pickup and from a DC component in the electron- 
ics that leaks into the data, we adopt a 'two-patch' ob- 
serving strategy that alternates between the target and 
a region of nominally blank sky. This strategy is sim- 
ilar to that used by CBI for the m easurement of CMB 
polarization (|Readhead et a/.l [20041 ). 

In this mode, a leading patch that contains our target 
is tracked for three minutes, and then a trailing patch of 
the presumably blank sky is tracked for the same period 
of time. The trailing patch is located at an RA 3"^ 10® 
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Fig. 3. — Modulation between the target (left) and the 
blank sky (center) successfully removes the strong contribu- 
tions from ground pickup and electronic offset. The difference 
map (right) clearly reveals the signal from cluster A2142. 

greater than the RA of the target. Since it takes 10 sec- 
onds for the telescope to slew between the fields, the two 
patches are observed over identical azimuth and eleva- 
tion tracks. During both tracks the platform polariza- 
tion angle is controlled so that each baseline corresponds 
to a fixed u-v mode. The recorded data during the two 
tracks should then contain the same contribution from 
ground pickup, which is removed by a subsequent differ- 
encing of the two tracks. This strategy requires that both 
the electronic offset and the ground emission are stable 
within 6 minutes, which is far sh orter as compare d with 
the measured time dependence (|Lin et al\ l2008f ) . The 
penalty for adopting such a two-patch strategy is the 
loss of efficiency by a factor of 4, where a factor of 2 
comes from doubling the observing time and the other 
factor of 2 from doubling the noise variance when dif- 
ferencing the two patches. We also note that there is a 
CMB component as well as possible combinations from 
Galactic foreground and extragalactic point sources in 
the differencing map though at a level much lower tha n 
the expected SZE of massive clusters ( Liu et a/]|2008| ). 
Nevertheless before observations we did check our fields 
for such contaminations. 

Figure [3] demonstrates the application of the two-patch 
observing strategy for SZE mapping. The left panel is 
an image constructed from the leading patch, contain- 
ing cluster Abell 2142 (see Sec. |4]for the formalism and 
procedure of data analysis). The middle panel shows the 
corresponding image from the trailing patch of blank sky. 
It is clear that both patches are contaminated by signals 
much stronger than the SZE (mainly from the ground 
pickup), and that the images are closely similar. The 
difference of the two patches reveals the cluster signal 
(right panel), whose signal strength is only about 1% of 
the contamination signal in individual patches. 

3.4. Calibration Events 

To calibrate these data, we observe planets roughly ev- 
ery 3 hours using the same 2-patch observing strategy. 
This time the two patches are each observed for four min- 
utes and separated by 4"^ 10® in RA. Such a separation in 
RA is equivalent to an angular separation of at least 57 
arcminutes, assuring that leakage of the planetary signal 
into the second patch is small, less than 0.1% of the plan- 
etary flux density according to the primary beam pattern 
measured in Sec. 13.11 The detailed formalism on how the 
data are converted to calibrated visibilities is presented 
in Sec. 14. 3[ and the det ermination of the planetary flux 
densities is discussed bvl Lin et al\ ([2008). Long-duration 
observations of the planets show that the system is sta- 
ble at night (between about 8pm and 8am) so that a 
planetary calibration every three hours controls the cal- 
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TABLE 1 

Observing log for AMiBA SZE clusters. 



Fig. 4. — Dirty images of Saturn (left) and the quasar 
J0325+2224 (right), both cahbrated by Jupiter using the for- 
mahsm presented in Sec. IH 

ibratio n error to be within 5% in gain and 0.1 radian in 
phase (|Lin et a/.i r2Q08). 

To test the overall performance of our observing, cal- 
ibration, and analysis strategy, we first used the Sun 
to successfully phase-calibrate an observation of Jupiter, 
and then Jupiter to phase and amplitude calibrate an 
observation of Saturn. The left plot of Figure [H shows 
a dirty image of Saturn calibrated by an observation of 
Jupiter made four hours later. The FWHM of this im- 
age is identical to that of the dirty b eam shown in F i^- 
ure [21 and the index of asymmetry (|Wu et al\ l2QQlf ) of 
the cleaned image is well below 5%. 

3.5. 2007 observations and SZE Targets 

After first light of the seven-element AMiBA in 
September 2006, which was a drift scan on Jupiter and 
whose image was quickly constructed, several months 
of start-up tests were used to study the instrumental 
properties, to fine tune our system, to study the ground 
pickup, and to find an optimal observing strategy. Dur- 
ing this period, Saturn, Mars, Venus, and the Crab Neb- 
ula were also observed for system tests. In February 
2007, when testing the two-patch observing strategy, we 
detected a first extragalactic object, quasar J0325+2224 
with a total (unresolved) fiux density of 550 mJy (right 
plot of Figure [4|) . 

After further fine tuning of the system, in April 2007 
we detected our first SZE cluster, Abell 2142 with a 
central SZE decrement of around -320 mJy per beam. 
During 2007 we observed a total of six S-Z clusters with 
0.091 <z< 0.322: A1689, A1995, A2142, A2163, A2261, 
and A2390. 

The observations were always at night time, when the 
system is stable. Due to the weather, system engineering, 
and observability of the targets, a total of roughly 150 
observing hours was spent in observing the clusters or 
their trailing fields over a span of 5 months, from April 
to August in 2007. During this period a lot of efforts 
were dedicated to the sy stem tests and tuning (|Lin et al\ 
120081 : iKoch et a/.ll2008a[ ). Table □ summarizes the point- 
ing directions, the observing and integration times, and 
the synthesized beam sizes, for each of the six clusters. 
Due to the two-patch observing strategy, which modu- 
lates between on- and off-source positions, the total time 
spent during cluster observations is double the on-cluster 
time listed in the table. The integration time on source 
(tint) is less than the actual observing time (tobs) because 
of the data trimming that will be discussed in Sec. 14.11 

The six clusters were chosen based on observational 
convenience, the expected flux density based on SZE ob- 



cluster 


pointing d 


r. (J2000) 


^obs 


^int 


syn. 




RA 


DEC 


(hr) 


(hr) 


res. 


A1689 


13^11.49™ 


-1°20.47' 


8.75 


7.11 


6.3' 


A1995 


14^52.84™ 


58°02.80' 


15.90 


5.56 


6.6' 


A2142 


15^58.34™ 


27°13.61' 


7.05 


5.18 


6.5' 


A2163 


16^15.57™ 


-6°07.43' 


7.80 


6.49 


6.6' 


A2261 


17^22.46™ 


32°07.62' 


19.15 


8.87 


6.4' 


A2390 


21^53.61™ 


17°41.71' 


16.05 


11.02 


6.4' 



Note. — Observing log for the AMiBA observations of six 
clusters of galaxies, including the pointing directions, the observing 
(tobs) cind integration (tint) times per baseline, and the synthesized 
resolution (syn. res.; the FWHM of the azimuthally averaged dirty 
beam). 



ser vat ions at other frequencies, and the existence of ex- 
tensive published X-ray data. The observational con- 
venience requires that the targets can be seen at the 
Mauna-Loa latitude of 19.54° N, at the night time when 
the system is stable (between about 8pm and Sam), and 
during the main observing period between April and 
August in 2007. Several cluster catalogs are compiled 
and studied for the expected SZE flux at AMiBA fre - 



[2002h . VSA (34 GHz; 



iGrego et al 120011: 


Reese et al\ 


Lancaster et all 


2005), and 



The redshift ran^e of z <0.3 is so chosen that massive 
clusters of a typical virial radius of few Mpc is resolvable 
by AMiBA, whose resolution {&) corresponds to a physi- 
cal scale of 1.6 Mpc at z = 0.3 in a flat ACDM cosmology 
with matter density parameter ft^ = 0.3, cosmological 
constant density parameter ft a = 0.7, and Hubble con- 
stant Ho = 70 kms~^Mpc~^. Among these six clusters, 
only A1689 and A2390 are known to have relaxed with 
a nearly isothermal profile. A1689 has a regular mor- 
phology in X-ray while A2390 is elongated. A2163, by 
contrast, is a merger with shocks and high-temperature 
regions, and A2142 is a merger with a strong frontal 
structure. 

Due to the missing-flux problem, most studies for the 
scientific implications of our results require the assis- 
tance of X-ray derived properties of the clusters. In 
Table [21 we summarize the critical X-ray derived pa- 
rameters for the AMiBA clusters, which will be used 
by o ur companion science papers (Huang e^a/. 20()8; 
iKoch et a/.ll2008bl : ILiu et a/.l 120081 : lUmetsu a/.ll2008l) 
The X-ray temperatures are from iReese et al (120021) 



for A2 163, A2261, A1689^and A1995. iMarkevitch et all 
(119981 ) for A2142, and lAllen I (|2000D for A2390. The 
power-law index P in t he isotherma l /3-mo del and the 
core radius Or are from IReese et al\ (120021) for A2163. 



A2261, A1689, and A1995jSanderson fc PonmanI (|2003l ) 
for A2142, and iBohringer all (|1998D for A2390. The 



X-ray core radii are mostly below one arcminute except 
for A2142 and A2163. These parameters are mostly X- 
emission weighted but not corrected for any cooling flows, 
which we further discuss in some science a nalyses such 
as the estimation of the Hubble constant Hq (jKoch et al\ 

4. ANALYSIS METHOD AND RESULTS 
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TABLE 2 

Parameters of AMiBA clusters derived from X-ray 
observations. 



TABLE 3 

Percentages of data flagged by each flagging step. 



Cluster 



A1689 
A1995 
A2142 
A2163 
A2261 
A2390 



0.183 
0.322 
0.091 
0.202 
0.224 
0.232 



(keV) 

9.7±lf 
8.82dzQ" 



u.ouy^Q QQ5 



0.74±^; 
0.674±o 
0.516±o 
0.6 



3 

81i 

m 

013 



(arcsec) 
26. 6±^--^ 
38.9±|-§ 
188.4±]p 

15.7±|| 
28.0 



We describe the formalism that we use to analyze 
AMiBA data from their raw format to the calibrated 
visibilities, images, and cluster profiles. AMiBA oper- 
ates with a bandwidth of 16 GHz centered at 94 GHz. 
Its four-lag correlators provide two pass-bands of 8 GHz 
processing dual linear polarizations. The challenge for 
the analysis of such data arises because the measured 
output of the instrument is not visibilities and because 
there are only two channels of wide bands causing con- 
siderable band-smearing effect in frequency. With the 
formalism presented here we successfully used the Sun 
to calibrate the image of Jupiter, and then used several 
planets (Jupiter, Saturn, Mars, Venus) to cross-calibrate. 
Jupiter was chosen as the primary calibrator for our sci- 
ence results. 

When we observe a general source field 5'(x) with a 
synthesized primary beam 5(x), the visibility along a 
baseline b of length is 



v''{f)=S{k)<E>B{k), 



(4) 



where k is the corresponding u-v mode, / = |k|c/(27rL^) 
is the observing frequency, a tilde denotes the Fourier 
transform of a quantity, and (g) denotes a convolution. 
The four lag-outputs of a correlator for the source field 
can thus be modeled as 

c'{rm) = ^y\Hf)R\f)x 

exp [2^(/ - fo)rmi + ^\f)i] df} , (5) 

where (m = 1-4) is the time delay for each of the four 
lags, R^{f) and 0^(/) are the instrumental frequency- 
dependent gain- and phase-responses respectively, /o is 
the frequency of the Local Oscillator (86 GHz), and /i 
and /2 indicate the upper and lower ends of the response 
frequency range i.e. [fij2] = [86, 102] GHz for AMiBA. 
Our analysis must invert Equation ([5]) to obtain the cal- 
ibrated visibilities v^{f). 

Figure [5] shows our first-light lag data c^(r^) of a drift 
scan on Jupiter. The fringing rate is determined by the 
angle between a baseline and the drift direction. 

4.1. Time-Domain Pre- Processing 

Initially the cluster data must be processed to remove 
baselines for which the system is malfunctioning. In ad- 
dition to checking the instrumental logs of receivers and 
correlators, nightly drift-scan data for Jupiter or Saturn 
are used to identify obvious problems. We then apply 
Kolmogorov-Smirnov (K-S) test to filter out datasets for 
which the noise appears non-Gaussian (|Nishioka et al\ 
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'2008). In the K-S test, a 5 % significance level is used 
throughout and shown to be efficient in identifying data 
sets with known hardware problems. Finally, we remove 
periods of data where occasional mount control problems 
appear. 

After this flagging, we difference the tracking data 
c^(t^; t) in the two-patch observation in order to remove 
ground pickup and the electronic DC offset. 4cr outliers 
in the temporal domain are clipped before the differenced 
tracking data are integrated over the 3-minute (4-minute 
for the calibration) tracking period to yield the lag data 
c^{Tm) = {c^{Tm]t))t' Thus cach two-patch observation 
yields a set of four lags c^(rm) {m = 1-4). We carefully 
verified that the noise of the lag data is white within the 
time scale of 10 minutes, assuring the scaling of the noise 
level with the two-patch integration time (|Nishioka etlill 
[2008h . 

Further flagging of 4(T outliers is applied in the visibil- 
ity space (see Sec. l4.3p . Table [3] shows the percentages of 
data that are flagged out by each step. On average the 
fraction of good data is about 60%. 

4.2. Lag-to- Visibility Transform 

We must invert Equation ([5]) to convert the four lags 
c^(Tm) to two visibilities v^{f). Ideally, if = r is con- 
tinuous, then we need only perform an inverse Fourier 
transform from the r domain to the / domain. How- 
ever r here is sampled discretely, and conservation of the 
degrees of freedom provided by four lag measurements 
implies that we can obtain at best two uncorrelated com- 
plex visibilities, with no redundancy. 

We perform the inversion first by dividing the 16- GHz 
frequency band into two, and [/(i,/2], each asso- 

ciated with one 'band- visibility', and define a 'visibility 
vector' as 

nv'iihJd])} 

5RK([/d,/2])} 

^{vH{h,f2])} 



(6) 



If we also write the four lags in column- vector form = 
c^i'^m) (^ = 1-4), then the visibility-to-lag transform 
resulted from Equation ([5j) is simply = U^v^, where 
is a 4 X 4 matrix. The four columns in can be 
easily obtained from Equation ([5j) as the four vectors 
when setting each component in to unity in turn 
while zeroing the other three components. Subsequently 
the inversion is 



= [U^] = VV. (7) 

Once given any lag data c^, we can use Equation ([7j) to 
construct the band- visibilities v^{[fi, fd\) and v^{[fd^ /2])- 
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Fig. 5. — First- light lag data c^(rm) of a drift scan on Jupiter in September 2006. Each panel represents one baseline, with an 
'Xy label at the bottom-right corner indicating the correlation between the antennae X and Y. X = 1 represents the central 
antenna, with the other six close-packed around it. In each panel, the four black curves and four red curves correspond to the 
lag outputs of the LL and RR cross-polarization modes respectively (see Sec. [2]). Here a low-pass filtering at 0.05 Hz is used to 
remove high-frequency noise. The horizontal time axis is the offset relative to the transit time. 



In principle, if R^{f) and ^^(/) are accurately known, 
then can be accurately determined so that there is no 
need for further calibration. However, in reality R^{f) 
and (/)^(/) are time-dependent so that is also time- 
dependent and thus needs to be calibrated for every ob- 
servation. Because the transform ([7]) is linear, any rea- 
sonable choice of fd^ R^{f) and should yield the 
same after calibration. Hence we make the simplest 
choice fd = (/i + /2)/2 = 94 GHz, R\f) = 1, and 
^^(/) = and then calibrate the data using the for- 
malism described in Section 14.31 This choice makes the 
inversion ([7j) exactly a discrete inverse Fourier transform 
that accounts for the band-smearing effect. An alterna- 
tive, and equally natural, choice might be to have equal 
power in the two frequency bands, with fd so chosen that 

jj^^ R\f)df = jjl R\f)df, but this is difficult to imple- 
ment since R^{f) is time-dependent. 

4.3. Calibration 

For calibration, we observe planets using the same 
2-patch observing strategy (Sec. 13. 4p . the same time- 
domain pre-processing (Sec. 14. ip . and the same lag-to- 
visibility transform (Eq. ([7j), Sec. 14. 2p . The underlying 
band- visibilities of these calibrators have the form 



V -I- 



Ml 


M2 




(8) 



where Mi and M2 are real constants, corresponding 
to the fluxes at different frequencies. However, the 
band- visibilities constructed from observational data us- 
ing Equation ([7j), where we have chosen R^{f) = 1 and 
/(/) = 0, will take the form 
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Since the involved operations are linear, there must exist 
a calibration matrix that corrects this discrepancy, 
i.e. 

' ' ' ^ 6(obs) 



6(thy) 



(10) 



We may model as 
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(11) 

where the parameters a\ and account for the gain- 
and phase-corrections relevant to the effects from R^{f) 
and (j)^{f) in Equation ([5]) respectively. By comparing 
the moduli and phases of vl^^^^"^ and vl^^^^"^ , we obtain 

^6(obs) 

(12) 



il = Ml /Ml' 



= M2/M^^°^'\ 



Finally, for a set of observed lag data c^, the calibrated 
band- visibility vector is 

= C^V^c^ (13) 

We emphasize that this is a linear transform, so different 
initial choices for R^{f) and (j)^{f) that generate differ- 
ent should lead to different but the same final 
calibrated band- visibilities v^. For each two-patch ob- 
servation, we obtain one for each cross-polarization 
mode and each baseline. For each baseline and each day, 
4<j outliers in |v^| are flagged. 



4.4. Noise Estimation 

In further analyses such as image making and model 
fitting for the cluster profile, a reliable estimate for the 
error in the band- visibilities derived above will be needed. 
First the 2-patch differenced lag data c^{t) = c^{rm]t) 
are modeled as the linear sum of the signal and noise, 
(9) c^(t) = s^{t) + n^(t). Since such data for clusters are 
dominated by e lectronic noise (i.e. n^ (t) s^(t)) and the 
noise is white (INishioka et a/.ll2008f ), the noise variance 
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TABLE 4 
Properties of AMiBA SZE Images 



cluster 


dirty 

size 


images 

flux 
(mJy) 


size 


cleaned 

flux 
(mJy) 


images 

S/N 
ratio 


scale 
(Mpc) 


A1689 
A1995 
A2142 


(6.10 
(6.80 
7.7' 


-217 
-167 
-320 


(5.70 
(6.80 
9.0' 


-168 
-161 
-316 


6.0 
6.4 
13.7 


(1.05) 
(1.91) 
0.92 


A2163 


7.8' 


-347 


11.2' 


-346 


11.7 


2.24 


A2261 
A2390 


(6.20 
7.4' 


-115 
-180 


(5.8') 
8.0' 


-90 
-158 


5.2 
6.6 


(1.25) 
1.78 



Note. — Basic properties measured from the dirty and cleaned 
SZE images, including the angular sizes (azimuthally averaged 
FWHM), the peak flux, and the S/N ratios. For cleaned images, 
the peak flux is corrected for the attenuation by the primary beam 
due to an offset to the pointing center. The last column 'scale' 
indicates the physical scale at z (see Table O corresponding to 
the angular size shown in the fourth column. Here we have as- 
sumed a flat ACDM cosmology with = 0.3, 17 a = 0.7, and 
= 70 km s~-'^Mpc~^. The brackets indicate that the cluster 
appears unresolved in the SZE image. 

associated with each time-integrated = {c^{t))t can be 
directly estimated as 

KY^([na')^((n^(t))?)-^i^^, (14) 

itt 

where #t is the number of temporal data in the time in- 
tegration. Then (crj)^ constitute the diagonal elements 

of the lag-lag noise correlation matrix = (n^n^ ). 
The off-diagonal elements are found to be less than 10% 
of the diagonal elements (Nishioka et al. 2008) so we ap- 
proximate them as zeros. 

In a similar fashion, the band-visibility obtained from 
Equation ([13]) can be modeled as a linear sum of the 
signal and noise, v^(t) = s^(t) + n^(^). Finally the noise 
correlation matrix for the calibrated band- visibilities is 

= (n>^^) ^ C^V^N^V^^C^^. (15) 

4.5. Image Making and Cleaning 

We construct the 'dirty' image of a cluster by making 
a continuous inverse Fourier transform of all the band- 
visibilities from all days and all baselines, with noise 
weighting based on the noise variance. The instanta- 
neous u-v coverage is improved in AMiBA operation by 
rotating the platform to eight polarization angles, giving 
a uniform angular interval of 7°. 5 between sampled u-v 
modes k (see Sec. 13. 2p . Figure [6] shows the dirty images 
of these clusters. The flux decrement at the center of 
each target is evident, as expected for the SZE signal 
at the AMiBA center frequency of 94 GHz (see Sec. [2]). 
The left part of Table [H summarizes the angular sizes (az- 
imuthally averaged FWHM) and the peak fluxes directly 
measured from these images. 

A dirty beam for each cluster data set is constructed 
using the same method, and the dirty image and 
dirty beam are then processed by a CLEAN proce dure 
(^Hogbom.1974) in MIRIAD (MIRIAD-ATNF 11200 8) to 
yield a cleaned image. Figure [71 shows the cleaned SZE 
images of the six AMiBA clusters, where the cleaned re- 
gions are indicated by the white circles, which are the 
FWHM contour of the primary beam. This process sig- 
nificantly reduces the convolution effects from the finite 
u-v coverage. The basic properties measured from the 



cleaned images are summarized in the right part of Ta- 
ble [H The apparent angular sizes, as compared with the 
synthesized resolution of about 6^ indicate that we have 
partially resolved clusters A2142, A2163, and A2390, 
while A1689, A1995, and A2261 appear unresolved in 
these images. The SZE signals observed here are at the 
level of few hundred mJy, indeed consistent with the the- 
oretical expectation A/sze ~ —200 mJy as discussed in 
Sec. [2 

In Table [H the S/N ratios are computed as the peak 
flux of the cleaned model (uncorrected for the primary- 
beam attenuation) divided by the RMS of the noise resid- 
ual map. We note that the noise level in these results 
are dominated by the instrumental noise, leading to the 
fact that the scaling among the peak flux, S/N ratios, 
and integration time is consistent with the noise equiv- 
alent flux (point-source sensitivity) of 63 mJyA/hi\ for 
the on-sou rce integration time in a 2-patch observation 
(|Lin et a/.| l2008). i.e. for a 2-hour observation, one hour 
per patch, the noise RMS is 63 mJy. Simulations also 
show that our analysis method (described in Sec. l4.21[0|) 
does not bias the SZE amplitude but induces a statisti- 
cal error at the few-per cent level. The calibration error is 
controlled within 5 % (|Lin et a/.ll2008f ). The systematic 
effect from CMB anisotropy is estimated to be at a sim- 
ilar level as the instrumental noise, and the point-source 
contamination causes a n underestimate of the SZE am- 
plitude by about 10 % (|Liu et a/.ll2008l ). 

4.6. Estimation of Cluster Profiles 

It is useful for science purposes to estimate parame- 
ters describing the structures of the SZE clusters from 
the visibilities. We employ a maximum-likelihood anal- 
ysis to estimate the band-visibilities taking the priors 
that all the visibilities represent the same circularly- 
symmetric source and that the visibility moduli within 
a frequency band are constant. Thus there are only six 
non-redundant band- visibilities {vp] p = 1-6} of six mul- 
tipole bands corresponding to the two frequency bands 
[fii f2] = 2) of three baselines lengths (see Sec. 13.21 
and Table [5]), i.e. the observed band- visibilities are mod- 
eled as 

r ^^'^^°"^^^^k(16) 

where xq is the SZE cluster center, xi is the pointing 
direction, 5(x) is the primary beam, and the integration 
accounts for the frequency band-smearing effect. These 
{vp} are equivalent to the band- visibilities at xi = xq 
and therefore are real numbers. They indicate the cluster 
profile in the visibility space. 

Markov Chain Monte-Carlo (MCMC) approach with 
Metropolis-Hastings sampling is used for this eight- 
dimensional likelihood analysis. For each cluster three 
MCMC chains of 200,000 samples are used. The results 
are given in Table [5] and Figure [HI For all clusters the ex- 
pected SZE decrement in flux density is evident. We note 
that although the primary CMB makes a no n- negligibl e 
contribution to our observed visibilities (^ Liu et a/.ll2008f ) , 
it appears with random phase in u-v space and thus will 
not bias our analysis here but rather contribute as part 
of the error bars. 
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Fig. 6. — Dirty images of the six SZE clusters observed by AMiBA. The white circles indicate the half- maximum contour of 
the primary beam (the field of view; 23^), while the blue regions at the bottom-right corners indicate the half- maximum regions 
of the noise- weighted synthesized beams 6'). The central decrements in flux density within the field of view provide strong 
evidence for the detections of SZE clusters. 



TABLE 5 

SZE Centers and Visibility Profiles for AMiBA Clusters 



cluster 


SZE center 


xo (J2000) 


Vi 


V2 


V3 


V4 


V5 


V6 




RA 


DEC 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


A1689 


13^11.41 ± 0.03"^ 


-1°20.7 ± 0.4' 


-123 ± 92 


-94 ± 89 


-566 ± 109 


-130 ± 108 


-414 ± 185 


85 ± 212 


A1995 


14^53.12 ± 0.06"^ 


58°02.6 ± 1.0' 


-234 ± 86 


-106 ± 80 


-241 ± 138 


-207 ± 126 


-53 ± 216 


-32 ± 200 


A2142 


15'"58.30 ± 0.03"^ 


27°13.8 ± 0.5' 


-508 ± 78 


-366 ± 76 


-140 ± 113 


-205 ± 117 


-9 ± 139 


-206 ± 162 


A2163 


16'^15.73 ± 0.04"^ 


-6°09.5 ± 0.5' 


-652 ± 110 


-332 ± 113 


93 ± 167 


-236 ± 158 


-457 ± 298 


-31 ± 260 


A2261 


17'^22.46 ± 0.06"^ 


32°08.9 ± 0.5' 


-34 ± 68 


-108 ± 67 


-160 ± 94 


-78 ± 99 


-412 ± 138 


-37 ± 163 


A2390 


21^53.72 ± 0.05"^ 


17°38.8 ± 0.7' 


-148 ± 87 


-185 ± 75 


-260 ± 107 


-273 ± 112 


-96 ± 141 


-44 ± 166 



Note. — The maximum-likelihood results for the coordinates of SZE cluster centers xq and band-visibility profiles {vp}. The values vi-ve 
correspond to the multipole bands of ^ = [1092,1193], [1193,1295], [1891,2067], [2067,2243], [2184,2387], [2387,2590] respectively (see Sec. [Ql 
Sec. 14.61 ). vi, vs, correspond to the frequency band of [86,94] GHz, and V2, V4, vq correspond to [94, 102] GHz. The three baseline lengths 
are 0.606 m (for vi, V2), 1.05 m (for V3, V4), and 1.212 m (for V5, ve)- 
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Fig. 7. — Cleaned SZE images of clusters observed by seven- element AMiBA. The white circles indicate the half-maximum 
contour of the primary beam (the field of view; 23^), and the blue patches at the bottom-right corners show the half-maximum 
regions of the dirty beams 6^). 



10 



Wu et al 




200 



-200 
-400 
-600 



1000 2000 
Multipole / 



1000 2000 
Multipole / 



Fig. 8. Visibility profiles {vp} of AMiBA SZE clusters. 
The data for the crosses are summarized in Table [5l while 
the boxes are the noise weighted average of the two frequency 
bands at each baseline length. 

We emphasize that these results for {vp} and xq are 
independent of cluster model. They can be further com- 
pared or combined with other experimental results, or fit- 
ted with a specific cluster model, for f urther study. Th e 
approach here is different from that in Liu et al} (|2QQ8l ). 
where an isothermal /3-model is directly fitted with the 
visibilities to estimate the central brightness /q, with the 
core radius Oc and the power index /3 taken from X-ray 
analysis. 

5. TESTS FOR SYSTEMATIC ERROR 

To verify that our detections of SZE clusters are real 
rather than from the instrument or from the foreground, 
we implement several tests. 

5.1. Differencing Maps 

A commonly used powerful test is the so-called sum- 
and-difference test, where the temporal data are divided 
into two subsets of equal size and then processed sep- 
arately. Figure [9] shows an example for A2142, where 
the two half-data images (left and middle) show a clear 
signal while their difference (right) reveals no signal but 
noise. Their average is very close to the overall dirty im- 
age shown in Figure [6l We have verified that such feature 
for the existence of signal is independent of the scheme 
for dividing the data into two halves. All data of the 6 
clusters have passed this test. 
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Fig. 9. — Both dirty images (left and middle) of A2142 
analyzed from two equally divided sub datasets show a clear 
SZE signal of flux decrement, while their difference (right) 
indicates no recognizable signal but noise. 



We carried out several other tests for systematics. 
Long-duration (12-hour) observations of planets were 
used to check the stability of the system; an indepen- 
dent analysis path produced results for A2142 consistent 
to high accuracy with those presented here; and several 
two-patch blank-sky observat ions showe d no signal above 
the expected CMB confusion (|Lin et a/.| [2QQ8). The noise 
properties of the lag data were investigated in detail indi- 
cating no s ystematics and no significant non-Gaussianity 
([Nishioka et al. 2008). The asymmetry of the antenna 
beams (Wu et al. 2001) was shown to be negligible. The 
pointing error is sma ll enough that it h as negligible ef- 
fects on our results (|Koch et al\ l2008a[ ) , and the radio 
alignmen t was tuned to im prove efficiency and tested for 
stability (|Wu et a/.ll2008b[ ). 

6. SCIENTIFIC IMPLICATIONS 

The scientific implications of the AMiBA SZE results 
presented here are further studied by companion papers. 
The Hubble constant is estimated by first deriving the 
angular diameter distances of th e clusters, and foun d to 
be i^o = 50t}^t^^ kms-^Mpc"^ (|Koch et a/.ll2"008bh in a 
flat ACDM cosmology with = 0.3 and (^a = 0.7. We 
also estimate the i ntegrated Compton parameter I2500 
(|Huang et a/.ll2008l ). which is the Compton y parameter 
integrated out to the angular radius at which the mean 
overdensity of the cluster is equal to 2500 times the criti- 
cal density of the universe at that redshift. Table [6] sum- 
marizes the results, as compared with the results from 
observations at 30 G Hz (OVRO/BIMA, deduced from 
iMorandi et al\ 12008 '. relativistic correctio n considered) 
and 145 GHz (SuZIE II, deduced from iBenson et al\ 
|2004 relativistic correction considered). Our results 
are consistent with those from OVRO/BIMA except for 
A1995, but seem to be systematically lower than those 
from SuZIE II. We also investigate the scaling relations 
between the I2500 from AMiBA data and various X-ray 
derived properties such as the gas temperature Te, to- 
tal mass M2500 and luminosity L^^. The scaling powers 
of these relations are consistent with the predictions of 
the self-similar model (|Huang et a/.ll20Q8l ). In the above 
studies, due to the missing flux problem the spherical 
isothermal /3 model is used with the spectral index f5 
and the core radius estimated by the X-ray data (see Ta- 
ble [2|), and the normalization calibrated by th e AMiBA 
SZE visibilities analyzed here (' Liu et a/]|2008[ ). 

We also perform a joint analysis of our SZE data with 
the weak gravitational lensing data from Subaru obser- 
vations (Umetsu et al. 2008). For the four clusters of 
A1689, A2142, A2261, and A2390, the two data sets are 
found to be in great agreement in morphology. Quanti- 
tative analysis even yields an estimation for the baryonic 
fraction of /b(< ^200) = 0.133 ± 0.027 and concludes 
that when compared wi th the cosmic baryo nic fraction 
r^b/^m = 0.171 ± 0.009 (jPunklev aI]l2QQ^ . 22 ± 16 % 
of the baryons are m issing from the hot phase of clusters 
(|Umetsu et a/.ll2QQ8l ). The morphological agreement be- 
tween our data and the lensing data and the consistency 
of /b, i^o, ^2500 and scaling relations with literature are 
encouraging, indicating that AMiBA is a reliable CMB 
telescope. 



5.2. Other systematic tests 
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TABLE 6 

Comparison of AMiBA SZE Results with Other 
Observations. 



Cluster 


1^2500 (xlO--^^ sr) 
OVRO/BIMA AMiBA 


SuZIE II 


A1689 


2.17 ±0.14 


3.13 ±0.96 




A1995 


0.71 ±0.06 


1.60 ±0.36 


A2142 




14.67 ±2.62 




A2163 


5.53 ±0.41 


6.32 ±1.10 


O.OU_Q 7Q 


A2261 


1.51 ±0.18 


1.36 ±0.71 


4 4f.+i.70 


A2390 




1.69 ±0.67 


o aq+0.56 



Note. — The integrated Compton parameter 12500 measured 
by AMiBA (86-102 GHz; Huang et al. 2008) as compared wi th re- 
sults from OVRO/BIMA (30 GHz; deduced from Morandi eToD 
1 20081 , relativist ic correction consid ered) and SuZIE II (145 GHz; 
deduced from IBenson et a/.l 120041 . relativistic correction consid- 
ered). 



We successfully detected six SZE clusters with the 
seven-element AMiBA in its compact configuration. The 
analysis method and results presented here mark a mile- 
stone for the AMiBA project and provides the first suc- 
cessful results for SZE clusters in the 3-mm band. These 
results are consistent with the published results based 



on data from other observations, and provide comple- 
mentary information at the same time. A number of 
tests show that the system performs as expected from its 
design, and we anticipate that t he current expansion to 
a 13-element system (Ho et al 2008b'; 'Wu a/."2008a') 
will boost its capability for the study of CMB cosmology. 
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